Spin Triplet Excitations for a Valence Bond Solid on the Kagome Lattice 
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One of the most promising candidate ground states for the quantum antiferromagnetic Heisenberg model on 
the Kagome lattice is the valence bond solid (VBS) with a 36-site unit cell. We present a theory of triplet 
excitation spectra about this ground state using bond operator formalism. In particular we obtain dispersions 
of all 18 triplet modes in the reduced Brillouin zone. In the bond operator mean-field theory, it is found that a 
large number of triplet modes are non-dispersive. In particular, the lowest triplet excitation is non-dispersive and 
degenerate with a dispersive mode at the zone center. Away from the zone center, the lowest triplet is separated 
from two other flat modes by a small energy gap. Quantum fluctuations are considered by taking into account 
scattering processes of two triplets and their bound state formation, which leads to a downward renormalization 
of the lowest spin triplet gap. The dispersion of the lowest triplet excitation in the VBS state is compared with 
the dispersive lower bound of the triplet continuum expected in competing spin liquid phases. Implications to 
future neutron scattering experiments are discussed. 

PACS numbers: 74.20.Mn, 74.25.Dw 



I. INTRODUCTION 

The quantum antiferromagnetic Heisenberg model on the 
Kagome lattice is a quintessential example of frustrated quan- 
tum magnets and therefore has been a subject of intense re- 
search activitie a'-^-^''*'^'Vi^i^i'°'''''^ . The ground state of this 
model, however, has been highly controversial generating a 
number of competing proposals. Recent proposals include a 
gapless spin liquid with Dirac fermionic spinons^'"^, a gapped 
spin liquid with bosonic spinons^^, and the valence bond solid 
(VBS) with a 36-site unit cell structure^^ii^. 

The identification of the ground state may be important 
to explain a series of recent experiments on herbertsmithites 
ZnCu3(OH)6Cl2 '3, '4, i5i 16, where spin- 1/2 moments reside 
on the ideal Kagome-lattice structure and interact with each 
other antiferromagnetically. It was found that the material re- 
mains paramagnetic down to 50 mK while the Curie-Weiss 
temperature is about 300 K'^i'^^i'^i'^i'^ . While this discovery 
is consistent with nonmagnetic ground states, the precise na- 
ture of the ground state still needs to be clarified by future 
experiments. 

Recent theoretical studies reveal that the energy differ- 
ence between various competing ground states is extremely 
smal l^'^''"'"i'^ , which makes it quite difficult to theoretically 
determine the true ground state. At the same time, this implies 
that even a small perturbation to the ideal nearest-neighbor 
Heisenberg model on the Kagome lattice may result in com- 
pletely different ground states. Given that small perturbations 
such as the Dzyaloshinski-Moriya interaction and impurities 
are indeed likely to be present in herbertsmithites'^, the iden- 
tification of the true ground state in this material becomes a 
complex issue. In this context, it is crucial to make testable 
experimental predictions for various proposed ground states, 
which then can be used to determine the true nature of the 
nonmagnetic state discovered in herbertsmithites ' 3, 1 5. 1 6, i? _ 

In this paper, we focus on the valence bond solid ground 




FIG. 1 : (Color online) Dimer covering pattern for the valence bond 
solid phase with a 36-site unit cell. The unit cell with 36 sites (or 
18 dimers) is indicated by the enclosing dotted line. Thick solid bars 
denote singlet dimers. We use numbers between 1 and 9 to represent 
surrounding dimers (light-blue bars) while core dimers (dark-blue 
bars) are denoted by numbers between 10 and 18. Red circles indi- 
cate the locations of perfect hexagons. Note that perfect hexagons are 
made of three neighboring dimers, (1, 4, 7) or (3, 6, 9). The dimers 
2, 5, and 8 become bridges connecting perfect hexagons. 



state with a 36-site unit cell (shown in Fig. [U and compute 
spin-triplet excitation spectra by using bond opertor formal- 
ism. The triplet excitation spectra can be directly measured 
in future neutron scattering experiments when single crys- 
tals become available. The eighteen dimers in a 36-site unit 
cell can be categorized into two groups: "core" dimers (num- 
bered 10 — 18 in Fig. [T]i including a "pin-wheel" structure 
at the center and "surrounding" dimers (numbered 1 — 9 in 
Fig. [Hi including the honeycomb-lattice structure of "perfect 
hexagons". As shown later, triplet excitations in the core 
dimers are highly localized and have flat dispersions. On the 
other hand, triplet excitations in the surrounding dimers can 
develop a dispersion by hopping around perfect hexagons. 



2 



Triplet excitation spectra are obtained by using the fully 
self-consistent mean-field theory in the bond operator repre- 
sentation. As shown later (see Figl3]and Fig|6]l, the main fea- 
ture of the triplet excitation spectra is that a large number of 
triplet modes are non-dispersive. The lowest triplet excitation 
is non-dispersive and degenerate with a dispersive mode at the 
zone center. Away from the zone center, the lowest triplet is 
separated from two other flat modes by a small energy gap. It 
is expected that the triplet excitation energies at high symmet- 
ric points such as F and X may be useful for comparison to 
future neutron scattering experiments on single crystals. 

The mean-field energy of the ground state is — 0.427J 
which is not too bad for a naive mean-field theory when com- 
pared to —0.438 J obtained from exact diagonalization of a 
36-site cluster^. On the other hand, the lowest spin-triplet gap 
is rather high, namely 0.795 J, which is expected to be signif- 
icantly reduced once quantum fluctuations are properly taken 
into account. It is shown later that the formation of two-triplet 
bound state can renormalize the lowest spin-triplet gap from 
0.795Jto0.622J. While itis still larger than 0.164 J obtained 
from exact diagonalization of a 36-site cluster'^, it is certainly 
in the right direction. Possible origins for the discrepancy are 
discussed later 

The rest of the paper is organized as follows. In section HIl 
physical motivations for studying the VBS state with a 36- 
site unit cell are provided. In this section, the bond operator 
representation is explained. In section |III] the bond operator 
mean-field theory is developed. Results for the spin-triplet ex- 
citation spectra and quantum fluctuation effects are presented 
in section |IV] Finally, in section |V] we conclude by making 
a direct comparison between the lowest spin-triplet excitation 
obtained from our theory and those from the spin liquid the- 
ories. In the spin liquid states, triplet excitations form a con- 
tinuum. The lower bound of such continuum, or the threshold 
energy, is given by the convolution of spinon-antispinon ex- 
citation spectra. Detailed measurement of the triplet spectra 
would be a key to the identification of the true ground state. 



II. MOTIVATION FOR VALENCE BOND SOLID 

In this section we discuss physical motivations for the va- 
lence bond solid (VBS) state with a 36-site unit cell structure. 
The VBS state with a 36-site unit cell was initially proposed 
by Marston and Zeng^ who envisioned the Kagome lattice as 
a honeycomb-lattice arrangement of perfect hexagons com- 
posed of resonant spin singlets. 

Some time later, Nikolic and Senthil** provided an argu- 
ment for the general validity of the 36-site VBS state in the 
Kagome-lattice antiferromagnet by a duality mapping. The 
main idea of this work is that the original Kagome-lattice an- 
tiferromagnetic model can be mapped to the fully frustrated 
Ising model on the dual dice lattice with transverse fields. 
Using a reasonable assumption about the magnitude of the 
transverse fields relative to the Ising coupling parameter, it 
is shown that the honeycomb structure of perfect hexagons is 
likely to be the ground state. Meanwhile, in a series expan- 
sion study by Singh and Huse'', it has been shown that the 



ground state energy is minimized when perfect hexagons are 
connected to their neighbors through empty triangles sharing a 
singlet dimer bond which in turn implies the same honeycomb 
structure of perfect hexagons as discussed in earlier works. 

Here we provide an alternative viewpoint on the origin of 
the valence bond solid state. We begin by asking what state 
may be the most stable dimer covering configuration on the 
Kagome lattice. We argue that the most stable configuration 
is the dimer covering which maximizes the number of "topo- 
logically perpendicular" spin-singlet dimers. This argument is 
motivated by the exact ground state of the antiferromagnetic 
Heisenberg model on the Shastry-Sutherland lattice, where all 
dimers are mutually perpendicular to each other^". We use 
the word "topologically perpendicular" since two neighboring 
dimers have the identical spin exchange energy as those ob- 
tained in the Shastry-Sutherland lattice when a spin belonging 
to a given dimer shares a triangle with its neighboring dimer 

It can be shown that the pinwheel structure at the core of 
the unit cell in Fig. [T] is the configuration maximizing the 
number of "topologically perpendicular" dimers. The rest of 
the dimer covering falls into the honeycomb array of perfect 
hexagons. While the final conclusion is exactly the same as 
before, there is an additional advantage over the previous ar- 
guments. Our point of view naturally leads to the fact that 
core dimers are decoupled from surrounding dimers. To see 
this, it is convenient to use a mathematical formalism which 
uses the spin-singlet degree of freedom as a natural building 
block of the VBS phase. The bond operator representation is 
such a formalism^ii22i21. 

Let us consider two neighboring S = \ spins, S/j and S^. 
The Hilbert space is spanned by four states which can be de- 
composed into a singlet state, \s), and three triplet state, \tx), 
\ty) and \tz)- Then, singlet and triplet boson operators are 
introduced such that each of the above states can be created 
from the vacuum |0) as follows: 
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To eliminate unphysical states from the enlarged Hilbert 
space, the following constraint needs to be imposed on the 
bond-particle Hilbert space; 



-tU - 



1, 



(2) 



where a — x,y, and z, and we adopt the summation conven- 
tion for the repeated index hereafter unless mentioned other- 
wise. 

Constrained by this equation, the exact expressions for the 
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spin operators can be written in terms of the bond operators. 
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(3) 



where Safij is the third rank antisymmetric tensor with Sxyz — 
1. 

When neighboring dimers are "topologically perpendicu- 
lar", spin-singlet contributions from both of the constituent 
spins of the neighboring dimer exactly cancel out. This can- 
cellation can be seen in the bond operator representation of 
spin operators in Eq. Q where it is shown that a pair of the 
spin operators within the same dimer have the opposite sign in 
the parts containing spin singlets. This results in the Hamilto- 
nian with no dispersive quadratic part for core triplets, which 
eventually gives rise to nine-fold degenerate flat bands plotted 
in Fig. |3]and|6]in Sec.|IV] This fact is completely consistent 
with conclusions from the series expansion studie a"i^^ . 

In conclusion, the honeycomb lattice of perfect hexagons is 
a stable dimer covering maximizing the number of "topologi- 
cally perpendicular" dimers, which in turn naturally suggests 
decoupling of core dimers from surrounding ones. Since the 
bond operator representation is a perfect theoretical frame- 
work for such situation, we use it to analyze the Kagome- 
lattice antiferromagnetic Heisenberg model. 

Since we expect qualitatively different behaviors between 
core and surrounding dimers, we introduce a different set of 
parameters for their singlet condensate densities, (si.„) = s„ 
and chemical potentials, /ii „. Here i denotes the position of 
the unit cell located at i and n indicates the dimer index inside 
the unit cell. Furthermore, it is expected that the dynamics is 
different for the dimers within perfect hexagons (the dimers 1 , 
3, 4, 6, 7, and 9 in Fig.[Tll and those bridging perfect hexagons 
(the dimers 2, 5, 8 in Fig.[Tll. Therefore, we introduce sc, He 
for nine core dimers, sh, IJ-h for six surrounding dimers of 
perfect hexagons, and ssr, fJ'Br for three bridging dimers con- 
necting perfect hexagons. Technical details of the bond oper- 
ator analysis are provided in the next section. Readers who 
are only interested in the results may directly go to Sec. |IV] 



where (i, n) denotes the position of the n-th dimer in the unit 
cell located at i. The chemical potential, /ii.„, is set to be /ic 
for core dimers, /i// for perfect hexagon dimers, and fiBr for 
bridge dimers. Similarly, the spin-singlet condensate density, 
s„, is set to be sc for core dimers, sh for perfect hexagon 
dimers, and SBr for bridge dimers. 

The total Hamiltonian, H = Hj + H^, can be written as 
follows: 

H ~ Neo + i?Quad, Coi-e + -ffQuad, SuiTounding + -ffQuartic, (6) 

where i?Quad. Core dcnotcs the quadratic part of the Hamilto- 
nian for core triplets while iJQuad, SuiTounding represents for sur- 
rounding triplets. In the above, N is the number of unit cells 
and 



£0 = 9 
+ 6 

+ 3 



Mc(l 
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- Sh) - - Jsjj 



-Js 



Br 



(7) 



It is important to note that the Hamiltonian does not have 
the coupling between core and surrounding triplets at the 
quadratic level. As shown in Sec. |IV] the absence of quadratic 
coupling is crucial for a complete decoupling between core 
and surrounding triplets. 

The quadratic Hamiltonian for core dimers is given by 



/J \ 

ffQuad.Core = ( T " A^C j ^1 XI (k)ina (k) , (8) 

^ ^ k n=10 

where a = x,y, and z. Also, the quadratic Hamiltonian for 
surrounding dimers is given by 



III. HAMILTONIAN AND THE MEAN FIELD THEORY 

We consider the following Hamiltonian for the Heisenberg 
model: 

i?j= 51 •^"■,r'S(r).S(r') (4) 

(r,r') 

where r indicates the original coordinate of the lattice, J^y = 
J within dimers, and Jr.r' ~ A J between neighboring sites 
belonging to different dimers. Utilizing the bond operator rep- 
resentation of spin operators, the Hamiltonian can be rewrit- 
ten solely in terms of bond particle operators. At this point 
the hard-core constraint among bond particle operators is im- 
posed via the Lagrange multiplier method; 

= ~ X! + ^InJhria ~ 1), (5) 

i,n 



^Quad, Surrounding — -f^Quad.O + -f^Quad.A + -^Quad.B + ^^Quad.C 

(9) 

where i?Quad.o denotes the quadratic Hamiltonian of triplet op- 
erators within a unit cell. iJguad.A, -f^Quad.B, and -ffguad.c de- 
scribe quadratic coupling between a given unit cell and neigh- 
boring unit cells separated by displacement vectors, r^, r^, 
and rc, respectively. These displacement vectors are defined 
as follows: 

ta ~ 4\/3ai:, 

Tb = 2\^ax — (Say, 

rc = vb - va ^ ~2V3ax - 6ay, (10) 
where a is the distance between nearest neighbor spins. More 
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explicitly, we get 



explicitly, we obtain 



^ k negn 



+ 



k nSgBr 

E E{[4Jk)t„„(k)+H.c.] 

(m,n)eG5 k 

+ [*L(k)tL(-k)+H.c.]} (11) 



with gn = {1,4,7,3,6,9}, gsr = {2,5,8}, and Gs 
{(1, 2), (2, 3), (4, 5), (5, 6), (7, 8), (8, 9)}. Furthermore, 



-ffQuad,A = 

_s%XJ ^ ^{[e^>^r.4jk)W(k)+H.c.] 

{m,n)eGA,H k 

+ [e''^"--4a(k)4a(-k)+H.c.]} 
_Mj:AJ ^ ^{[e--4„(k)t„„(k)+H.c.] 

(m,n)eGA,Br k 

+ [e'''--iL(k)tL(-k)+H.c.]}, (12) 

where 

Ga,h = {(1,4),(9,6)}, Ga,b. = {(1,5),(9,5)}. (13) 

-ffQuad,B and -ffQuad,c Can be obtained from iJguad.A by re- 
placing (i) rA by tb and rc and (ii) (G^,// , GA.Sr) by 
{Gb,h, GB,Br) and (Gc,i?, Gc,Br), respectively. Here, 

Gb,h = {(9, 3), (7, 4)}, GB,Br = {(8, 3), (8, 4)}, 
Gc,H = {(7,l),(6,3)}, Gc,Br = {(7,2),(6,2)}. (14) 



The quartic interaction part of the Hamiltonian is analyzed 
in the Hartree-Fock mean-field theory via quadratic decou- 
pling. The resulting mean-field Hamiltonian is written as fol- 
lows: 



-f^Quartic — -ffQuarticO + -f^Quartic.A + -f^Quartic.B + -f^QuarticC 

(15) 



where, similar to the quadratic counterparts, -ffguartic.o denotes 
the quartic Hamiltonian of triplet operators within a unit cell. 
^^QuarticA, ifQuarticB, and i?Quartic,c describe quartic coupUng 
between a given unit cell and neighboring unit cells separated 
by displacement vectors, r^, r^, and rc, respectively. More 



-ffQuarticO — -^^JiQc ~ + QcS ~ ^Cs) 
+ 2^J{Q\ ~ Ph) + 2^J{Q%r1 ~ PbtI + Q\r2 ~ PBr2) 

+ T ^ E{[^^-i*-«('^)*-W+H-c.] 

{m,n)eGs k 

-[0Brl4a(k)4a(-k)+H.C.]} 

{m,n)eGc k 



-[Qc4a(k)tL(-k)+H. c.]} 
+ T ^ E{[^c:5iL(k)t„a(k)+H.c.] 

{m,n)eGcs k 

-[Qcs4a(k)tL(-k)+H.c.]} (16) 



where 

Gc = {(10, 11), (11, 12), (12, 13), (13, 14), (14, 15), 
(15, 10), (10, 16), (12, 17), (14, 18)}, 
Gcs = {(16, 1), (16, 9), (17, 4), (17, 3), (18, 7), (18, 6), 
(11, 2), (13, 5), (15, 8)}. (17) 

One also gets 



H, 



QuarticA 

AJ 

T 



E E {[J'He*-4«(k)W(k) + H. c] 

{m,n)eGA,H k 

-[QHe*--4a(k)iL(-k)+H.c.]} 
+ T ^ ^ {[PB.2e''^--tL(k)t„a(k) +H. c] 

{m,n)^GA,Br k 

- [QB.2e*--4a(k)iL(-k) +H. c.]}, 

(18) 

Similar to the quadratic case, i?Quartic,B and HQuartic.c can be 
obtained from -f/Quartic,A by replacing (i) by and rc and 

(ii) iGA,H,GA,Br) by {GB,H,GB,Br) and {Gc,H,Gc,Br), 

respectively. 

The above mean-field order parameters, Ph, Qh, Psri, 
Qsri, PBr2, Qsr^, Pc, Qc, Pcs, and Qcs, are determined 
by solving a coupled set of ten self-consistency equations. In 
other words. 



Pj — (^i na^j'"^*^) ' ^7 — (^i,na^j,ma) 



(19) 



where (i,j) indicates the positions of the neighboring unit 
cells and 7 e (G, H, Brl, Br2, CS), both of which are re- 
lated to the dimer index pair, {n,m). To be specific, if 
{n,m) G Gc, 7 = G andj = i. If {n,m) & Gs, 7 = Brl 
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FIG. 2: (Color online) Mean-field order parameters, P and Q, as 
a function of A. The top figure depicts dimer covering of the va- 
lence bond solid state in the vicinity of a perfect hexagon, which 
also shows how various Q's are defined. Note that each Q describes 
an off-diagonal correlation between two dimers connected by corre- 
sponding arrows. Diagonal correlations, Ph, Pbti, and Psr2, are 
defined in the same way. It is interesting to note that Qbti {Pbti) 
and QBr2 (PBri) have opposite signs. 



and j = i. On the other hand, if {n, m) e Ga,h, 7 = 
and j — i = ta- Those for Gb.h and Gc.h are defined sim- 
ilarly. Also, if {n,m) £ GA,Br, 7 = Br2 and j — i = ta- 
Those for GB,Br and Gc,Br are defined similarly. Finally, if 
(n, m) e Gcs, 7 — CS and j — i. 

In physical terms, Pc and Qc describe the diagonal and 
off-diagonal triplet correlations between core dimers, respec- 
tively. In the case of surrounding dimers, six order param- 
eters are introduced. First, Ph and Qh denote correlations 
between dimers within perfect hexagons. The other four or- 
der parameters, (PBrijQBri) and {PBr2,QBr2), represent 
correlations between a bridge dimer and two neighboring 
dimers belonging to the nearby perfect hexagons. In partic- 
ular, {PbtI, Q Bri ) describes the correlations between a given 
bridge dimer and the dimer lying inside the nearby perfect 
hexagon in its parallel direction. {PBr2,QBr2) denotes the 
other correlations. A schematic diagram showing the defini- 
tion of these order parameters is provided in the top panel of 
Fig|2] Finally, Pes and Qcs indicate the correlations be- 
tween core and surrounding dimers. 

In order to determine the ten mean-field order parameters 
defined above, one needs to compute yet another six unknown 
parameters which are three spin-singlet condensate densities. 



Sc. sh, and SBr, and three chemical potentials, /ic, /.ih, and 
liBr- The hard-core constraints for bond operators provide 
three equations: 

18 



n=10 k 



4. = i-3^ E J2(*LimM). 



(20) 



n£gBr k 



On the other hand, the energy minimization with respect to 
spin-singlet condensate densities generate the other three nec- 
essary equations: 



Bsc 



0, 



dSn 



dso 



ds 



H 



dsBr 



0. 



(21) 



Now the ground state energy and excitation spectra can be 
obtained as follows. For this purpose, the mean-field Hamil- 
tonian is rewritten in the following compact fashion: 



F = i V A^MA^ - -TrM, 



(22) 



where 



Aj.(k) = [tL(k), . . . , 4„(k), <i„(-k), . . . , ii8a(-k) 



(23) 



and the matrix, M, is determined by simply reorganizing the 
mean-field Hamiltonian. 

It is important to note that our Hamiltonian describes dy- 
namics of boson operators, in which case obtaining normal 
modes of the Hamiltonian is not equivalent to diagonalizing 
the matrix M in the above. Instead, we need to consider the 
following eigenvalue problem^"*: 



where 



I 
I 



(24) 



(25) 



with I being the 18 x 18 identity matrix and 



(26) 



with 77* = (771, ■ • ■ , ms.) and ^* = (Ci, ■ • . , Cis)- This differ- 
ence between fermionic and bosonic problems fundamentally 
originates from the difference in their operator commutation 
relations. 
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IV. RESULTS 



A. Ground state energy 



The ground state energy per unit cell is obtained by solving 
the eigenvalue equation described in the preceding section. In 
other words. 



£0 



^•HQc ^ -Pc + Qcs ~ ^cs) 



2 

3 
2 

3, .,^2 



XJ{Q% 



Pi) 



-\J{Q 



Brl 



3 1 



p2 

^Brl 



2N 



k n=l 



Qi 

18 / J 

y V4 



P'Br2) 



fJ-H 



(27) 



where cun {n — 1, . . . , 18) are triplet eigenenergies obtained 
from diagonalization. As mentioned previously, the ground 
state energy is minimized with respect to spin-singlet con- 
densate densities. Note that the minimization process is 
performed simultaneously satisfying the ten self-consistency 
conditions and three hard-core constraints. 

Now let us discuss numerical results. When quartic interac- 
tions are completely ignored, the ground state energy per site 
is found to be —0.414 J at A = 1. The inclusion of quartic in- 
teractions lowers the ground state energy to —0.427 J per site, 
which can be favorably compared with —0.438 J from exact 
diagonalization^. This result is quite encouraging given that 
it is obtained from a naive mean-field theory. This perhaps 
indicates the robustness of the VBS state. 



B. Triplet dispersions without quartic interactions 

We now discuss the energy dispersion of triplet excitations. 
In this section quartic interaction terms are ignored. The anal- 
ysis of quartic interaction effects is relegated to the next sec- 
tion. One of the reasons why we first focus on the quadratic 
part of the Hamiltonian is that the overall structure of the dis- 
persion is well captured even at the quadratic level. The inclu- 
sion of quartic interaction terms modifies the dispersion only 
quantitatively. 

As seen in Eq. ([8]), core triplets have no dispersions. In 
this situation the ground-state energy minimization condi- 
tion immediately leads to the conclusion that sc — I and 
He = —0.75 J, which implies that in the core part of the unit 
cell triplet fluctuations are entirely absent and triplet excita- 
tions are completely localized. As mentioned previously, the 
reason for this complete localization has to do with the dimer 
covering structure of the valence bond solid phase. The key 
fact is that every core dimer is "topologically perpendicular" 
to its neighboring dimers in the sense that each spin belonging 
to a core dimer shares a triangle with its neighboring dimer. 
This results in the Hamiltonian with no dispersive quadratic 



2.0 



(o(k)/J 



1.0 



K 



K 



FIG. 3: (Color online) Energy dispersion of the triplet excitations 
when quartic interaction effects are ignored. The dotted line de- 
scribes the nine-fold degenerate triplet states from core dimers. The 
solid lines denote the nine triplet states from surrounding dimers. 



part for core triplets, which gives rise to nine-fold degenerate 
flat bands plotted as a dotted line in Fig. [3] 

Now let us turn to the remaining nine triplets in the 
surrounding dimers. After diagonalizing the corresponding 
bosonic Hamiltonian via the method described in the preced- 
ing section, mean-field parameters, SH,SBr, I^h, and jiBr, are 
self-consistently determined. Fig. [3] shows the resulting dis- 
persion of nine triplet excitations at the isotropic exchange 
limit of A = 1. It is convenient at this point to classify char- 
acteristics of these eigenmodes using group theory. Since the 
unit cell possesses the C3 point group symmetry, the states at 
the r point can be decomposed with respect to irreducible rep- 
resentations of the C3 point group. A representation, R, can 
be decomposed in the basis containing nine triplets localized 
at each dimer as follows: 



(28) 



where A is an one-dimensional representation and _E is a two- 
dimensional one. (See Ref. IH for standard conventions and 
notations.) 

The lowest flat mode belonging to the E irreducible repre- 
sentation is degenerate with a dispersive mode at the F point 
{El mode in Fig|3]l. There are two other flat modes {Ai and 
A2) which are separated from the lowest flat mode by a small 
gap. These states have A characters and are completely local- 
ized within perfect hexagons. The fifth dispersive mode (A3) 
is another A mode. The remaining four modes belong to the 
E irreducible representation and make two doubly degenerate 
states at the F point {E2 and E^). Real space representations 
of the five flat mode eigenvectors are shown in Fig.|4]and|5] 
Details of these eigenvectors for non-dispersive modes are fur- 
ther analyzed in the next section where their explicit forms are 
also presented. In the next section we study quartic interaction 
effects via the self-consistent mean-field theory. 



C. Triplet dispersions witli quartic interactions 

In this section we solve the total Hamiltonian contain- 
ing contributions from all eighteen dimers in the unit cell. 
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(a) . (b) 
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FIG. 4: (Color online) Real space representation of the eigenvectors 
of the flat A mode at the F point. The size of ovals and the sign inside 
indicate relative weights and phases of the amplitude at each dimer, 
respectively, (a) One of the flat A modes with the energy eigenvalue 
SIa. Triplets around each perfect hexagon have the same phases and 
weights. Two different flat A modes are distinguished by relative 
phases between neighboring perfect hexagons, (b) Local geometry 
around a perfect hexagon. Here red dots indicate the position of R- 
spins within a dimer. See Eq. l[3} for the definition of the R-spin. 



(a) 




(b) 



4 5 ^6 
A E2 8\ 



(C) 




FIG. 5: (Color online) Real space representation of the eigenvectors 
of the flat E modes. To avoid degeneracy eigenvectors are computed 
slightly away from the F point to the direction of the K point, for (a) 
El mode, (b) E2 mode, and (c) £3 mode. Note that the size of ovals 
and the sign indicate relative weights and phases of the amplitude at 
each dimer, respectively. 



The fully self-consistent calculation, however, shows that 
core dimers are still completely decoupled from surrounding 
dimers and remain flat in energy. In other words. 



Pr 



P, 



cs 



Ics = 0, Sc = 1, 



-0.75J. 
(29) 



On the other hand, surrounding dimers experience band 
renormalization due to finite triplet correlations, P and Q. Or- 
der parameters for the triplet correlations are plotted as a func- 
tion of A in Fig|2l The order parameters are shown to have the 
largest value inside perfect hexagons which in turn implies 
that vacuum fluctuation is strongest at this location. Triplet 
correlations between bridge and perfect hexagon dimers are 
also quite large, about 86% of those between perfect hexagon 
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FIG. 6: (Color online) Energy dispersion of the triplet excitations 
when quartic interaction effects are included in the self-consistent 
mean-field theory. Away from the zone center, the lowest flat mode 
is separated from doubly-degenerate flat A modes by a small gap, 
Ao = 0.069 J. Eigenenergies of the flat modes are shown on the 
right hand side of the plot, where subscripts indicate the irreducible 
representation of each state at the F point. 



dimers. It is interesting to note that QbtI (PBri) and QBr2 
{PBr2) are the same in magnitude but are the opposite in 
sign. In Fig|6] we plot the triplet dispersions obtained from 
fully self-consistent parameters at A = 1. The overall struc- 
ture is identical to that obtained by ignoring quartic interaction 
terms. 

One of the most interesting characteristics of the triplet 
dispersion is the existence of a large number of flat bands. 
Among the nine modes coming from surrounding dimers, five 
states have no dispersion. (If one includes the flat modes orig- 
inating from core dimers, fourteen states out of eighteen are 
flat in energy.) These five flat modes are categorized in terms 
of the C3 point group symmetry at the T point. In the right 
hand side of Fig|6] each mode is indicated by , ,^Ei, 
r^Bj, and He^, respectively (these are also energy eigenval- 
ues). Note that ft a = = and fi^j < fi^j < ri^,. In 
the above notations subscripts denote irreducible representa- 
tions to which each flat mode belongs. 

It is interesting to note that the eigenvalue equation can be 
solved exactly for the flat modes. The precise analytic expres- 
sions for the flat mode eigenvalues are given as follows: 



ft El — fi-, flEs — 



(30) 



where 



al^ + [an + buf - qjj - 12((?L - ± ^ 

(31) 



- 24(g|, - 6|,) [a|, + {an + bnf - ql] 
+ 2AaBr {qsr + bsrTiaH + bH - qn) 
+ {qsr - bBr)^{aH + bn + qn) , 



(32) 
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and 



an 



J 

4 



4 



XJ 



AJ, 



IH ^ —{sh +Qh), QBr = —{sHSBr + QBr2)- 



(33) 



It is important to notice that the eigenvalues, Ha^ , , and 
fig^ , are completely determined by the parameters defined in- 
side perfect hexagons. In fact, the same energy eigenvalues 
are obtained by applying the bond operator theory to a sin- 
gle isolated perfect hexagon using the mean-field order pa- 
rameters, fiH, sh, Ph, and Qh- This means that for the flat 
A and E2 modes a simple product of perfect-hexagon eigen- 
states becomes that of the full Kagome lattice, at least in the 
self-consistent mean-field theory. 

As mentioned in the previous section, the flat A modes are 
completely localized within perfect hexagons, /. e. , there is 
no weight for the bridge dimers. Here triplets around each 
perfect hexagon have the same weight and phase while the 
Ai and A2 modes are distinguished by a relative phase differ- 
ence in neighboring perfect hexagons. The E2 mode is also 
localized within perfect hexagons. The difference between 
the flat A and E2 modes originates from the fact that they 
represent different eigenstates of the isolated perfect-hexagon. 
Schematic diagrams for the real space representation of these 
flat modes are provided in Fig.|4]and|5] which show the rela- 
tive weight and sign of the eigenvector amplitudes at each real 
space dimer location. 

To understand the flat nature of the A modes, one needs to 
examine the local geometry around a perfect hexagon which is 
plotted in Fig.|4jb). As seen in Eq. (O, two constituent spins, 
the R and L spins, within the same dimer are distinguished by 
the sign difference in the part containing spin singlet opera- 
tors. This is basically due to the odd parity of the spin singlet 
under the inversion with respect to the center of the dimer 
The i?-spins are denoted by red dots in Fig. HJb). Because 
of the three-fold rotational symmetry, every spin from bridge 
dimers is simultaneously connected to both R and L spins 
of the perfect hexagon dimers. Due to the above-mentioned 
sign difference, hopping amplitudes between perfect hexagon 
and bridge dimers cancel, leading to the complete localization 
within perfect hexagons. The flat E2 mode can be understood 
similarly. 

We now present precise analytic expressions for the eigen- 
vectors of all five flat modes in Table J] The prefactors, tm, 
in Table I] denote the relative magnitudes of the hole (parti- 
cle) component in the Bogoliubov quasiparticles (quasiholes), 
which are given as follows: 
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TABLE I: Eigenvectors of the five fiat modes. Among the 36 com- 
ponents of the full eigenvector, = (771, . . . , 7718,1^1, . . . , 
only the components for surrounding dimers are shown. Note that 
core dimers are decoupled from surrounding ones. In the above, 
21 — exp (— ik ■ ta), 22 = exp (— ik ■ rs), and (ri, r2) represents 
any pair of two complex numbers satisfying rir2 = —1. The prefac- 
tor, Tm, depends on the mode index, M £ Ai, A2, Ei, E2, and E3. 
Explicit expressions for tm as well as ci , C3 , di , and are provided 
in the text. 



qni^Ei ^ "-Br) + l'2aBrbBrqBr 



TE, 



qH{^%^ - a|^) + \2aBrbBrqBr 

- a%^){VLE3 - an - bn) 



qni^E^ - flsr) + l^GBrbBrqBr 
{(^^Esiblr -<&r)+ QaBriblr + <llr)} 



qni^ 



^Br) + I'ittBrbBrqBr 



(34) 



Note that ta-^, and depend only on the perfect 
hexagon parameters. The coefficients, ci (di) and C3 (d^) 
indicate the relative magnitude of coupling between perfect 
hexagon and bridge dimers in the Ei and E^ modes, respec- 
tively. More explicitly, one obtains 



Cl 



C3 



^Ei — ^Br 
bBr + TE^qBr 

VLe3 — a,Br 
bBr + TE^qBr 



di 



— Hei — a-Br 
bBr + qBr/TEi ' 

— i^Ba — a,Br 
bBr + qBr/TE3 



(35) 



Note that the eigenvector amplitudes in bridge dimers, /. e. , 
('72,?75,??8) and (6,^5,^8), are zero for the flat Ai, A2 and 
E2 modes, which confirms that they are completely localized 
within perfect hexagons. 

While the amplitudes in bridge dimers do not vanish, the 
nature of the Ei and i?3 modes is rather similar to that of the 



E2 mode. As one can see from Table U the amplitudes in- 
side perfect hexagons for the Ei and modes are precisely 
identical to those for the E2 mode. This identity is funda- 
mentally due to (i) the odd parity of the singlet operator and 
(ii) the three-fold rotational symmetry of the lattice. Coupling 
to bridge dimers splits triple degeneracy by lowering the en- 
ergy for the El mode and increasing it for the E^ mode. The 
difference from the E2 mode case is that hopping amplitudes 
between perfect hexagon and bridge dimers do not cancel, but 
instead they add. 

It is interesting to note that our Ei mode is fully consis- 
tent with the lowest triplet excitation obtained in the series 
expansion study by Singh and Huse'^. In fact, after properly 
redefining the unit cell convention, it can be shown that our 
El mode becomes precisely equivalent to the lowest energy 
eigenstate obtained in first order of A under the conditions that 
(i) all off-diagonal coupling terms are ignored, (ii) all quartic 
interaction terms are ignored, and (iii) perfect hexagon and 
bridge dimers are physically identical, i. e., sh = ssr and 
f^H = fJ-Br- Considering that our Ei mode does not change 
abruptly upon relaxing these conditions or restoring the full 
self-consistency, one may expect that our Ei mode is indeed 
adiabatically connected to the lowest energy eigenvector ob- 
tained in the previous series expansion study and perhaps 
the corresponding eigenvector in yet-to-be-studied higher or- 
der series expansion. 



D. Quantum fluctuations 

Up to now, triplet interactions are treated within the self- 
consistent mean-field theory. To investigate effects of quan- 
tum fluctuations, we need to go beyond the mean field theory. 
There are three classes of order parameters in the mean-field 
theory, which can be affected by quantum fluctuations. These 
are chemical potentials, /i, diagonal correlation parameters, 
P, and off-diagonal correlation parameters, Q. Spin-singlet 
condensate densities, s, which are the remaining variational 
parameters, can be determined by minimizing the ground state 
energy after quantum fluctuation effects are incorporated for 
the above-mentioned order parameters. Below we check how 
each of these order parameters is affected by quantum fluctu- 
ations. 

The chemical potential is the Lagrange multiplier for the 
hard-core constraint. In the usual weak-coupling limit where 
the singlet nature of the ground state is robust, the dominant 
contribution comes from the hard-core constraint, Eq. (|2]l, 
which can be conveniently implemented by an infinite on-site 
repulsion between triplets-^. When the triplet density is low, 
this hard-core constraint can be treated by summing ladder 
diagrams for the scattering vertex. The complete ladder dia- 
gram summation for the full lattice is rather complicated. For- 
tunately, however, in our system the most important aspect of 
the lowest excitation is determined by the nature of eigenstates 
inside a single perfect hexagon. Therefore, we expect that the 
dynamics inside a single perfect hexagon is a good indicator 
for the full lattice as far as the lowest excitation is concerned. 
It is shown that there is little difference between the chemical 




FIG. 7: Renormalization of the triplet pair emission (or absorption) 
amplitude by quantum fluctuations. Vertices come from the singlet 
5 = scattering channels of two triplets. 



potential obtained from the mean-field theory and that from 
the ladder diagram summation using a single perfect hexagon. 
For this reason, we ignore the effects of quantum fluctuations 
on the chemical potential. 

Next, we consider the effects of quantum fluctuations on the 
diagonal correlation order parameters, P, and the off-diagonal 
correlation order parameters, Q. From our mean-field analy- 
sis, it is shown that the shear size of the diagonal correlation 
parameters is almost one order of magnitude smaller than that 
of the off-diagonal counterparts (See Fig. |2|i. It is thus ex- 
pected that, while P is to be certainly renormalized by quan- 
tum fluctuations, the overall size of its renormalization should 
be much smaller than that of Q. Thus we focus on Q below 
while ignoring fluctuation effects on P. 

Off-diagonal correlations are enhanced when there are 
strong fluctuations toward the formation of two-triplet bound 
states, which is caused by the attractive interaction between 
nearest neighbor triplets^^'^^. Effects of these fluctuations on 
the singlet/triplet spectrum can be captured by considering 
successive particle-particle scattering processes which renor- 
malize the triplet pair emission/absorption amplitudes, /. e. , 
the coefficients of tH"^ + H. c. terms-''. To be concrete we 
describe below how the corrections in these coefficients are 
computed. 

In the self-consistent mean field calculation as described in 
the SecHn] the bare pair emission/absoprtion amplitude, 
is renormalized to B^^^ where 



^ I SH 



shsbt 



if7 = BrlorBr2, 



(36) 



and i?*^^ — B^ + Q^. To go beyond the mean-field theory, 
we consider scattering of two triplets: a + f3 ^ fj. + where 
triplet spin indices, a, /3, /i, and ly, belong to {x, y, z}. From 
the quartic interaction terms, one can get the following bare 
scattering amplitude: 



(37) 



which shows that in the singlet 5 = channel the scattering 
amplitude is given by = ^SapSf^^Vap^f^ty = 

We now evaluate the ladder series for scattering processes 



in Fig.|7] which renormalize B^^^ as follows: 



B. 



MF 



B^ — 



MF 



+ AB^ 



(38) 
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where 



and 



AD _ _ 7 2_ 



(39) 



"7 = E/ (k,^)Gi^^(-k,-^). (40) 

k 

In the above G^lf and G^m ^re the mean-field Green's func- 
tion for triplets in the n- and m-th dimer location in the unit 
cell. Conventions for the relationship between 7 and {n, m) 
are the same as those for P and Q in Eq. ( fT9l ). 

The inclusion of the above quantum fluctuation corrections 
modifies the triplet Hamiltonian matrix which, after diago- 
nalization, leads to a reduction of the lowest spin gap from 
0.795 J to 0.622 J. Our prediction for the lowest spin gap is 
still larger than 0.164J obtained from exact diagonalization 
of a 36-site cluster or 0.08 ± 0.02 J from a recent 7th order 
series expansion result (Note that the value from series expan- 
sion becomes 0.2 J in the same finite 36-site cluster)i2,. This 
may suggest that quantum fluctuations beyond what we have 
considered may be necessary to reach quantitative agreement. 
At the same time, finite-size effects also need to be examined 
very carefully. 



V. DISCUSSION 



Kagome-lattice system shows magnetic ordering when k is 

larger than 0.53^^. Results for the lower bound of the two- 

spinon continuum are plotted in Fig. [8] 
1.0 I 




kx 



FIG. 8: (Color online) Spin- triplet excitation spectra for three can- 
didate ground states of the Kagome-lattice Heisenberg antiferromag- 
net: (i) the valence bond solid (VBS) with a 36-site unit cell, (ii) 
the U(l) Dirac spin liquid and (iii) the Z2 spin liquid. In the case 
of the VBS state we plot the renormalized spectrum of the lowest 
triplet mode including quantum fluctuation effects as described in 
Sec. lIVDl For these two spin liquid states we plot the lower bound 
of the two-spinon continuum. The first Brillouin zone for the spin 
liquid states which has a 3-site triangular unit cell is depicted on the 
right hand side. The small dotted hexagon shows the first Brillouin 
zone for the VBS state with a 36-site unit cell. Two red dots lying 
in the y-axis indicate the positions (gauge-dependent) of the Dirac 
nodes for the U(l) Dirac spin liquid state. 



The Kagome-lattice antifeiTomagnetic Heisenberg model 
has been generally regarded as one of the most geometrically 
frustrated spin systems in two dimension. Because of this, 
the Kagome lattice has been considered as a promising can- 
didate system for realizing exotic quantum spin liquid ground 
states. In particular, there are two spin liquid states that have 
received much attention lately: (i) the U(l) Dirac spin liquid 
state suggested in a projected wave function study by Ran et 
al— and (ii) the Z2 spin liquid state obtained in a bosonic 
large-A^ Sp(A^) approach by Sachdev-'. 

To explicitly compare the lowest spin-triplet excitation en- 
ergy of the VBS state with those of the above-mentioned spin 
liquid states, we compute the lower bound of the two-spinon 
continuum in spin liquid phases. The lower bound (or the bot- 
tom of the continuum) is given by 

w(k) = Minp [e(k/2 + p/2) + e(k/2 - p/2)] , (41) 

where e(q) is the single spinon excitation energy with mo- 
mentum q. An assumption behind this calculation is that 
spinons themselves are weakly interacting. 

For the U(l) Dirac spin liquid state we use the tight-binding 
spinon Hamiltonian with the mean-field order parameter for 
spinon hopping x — 0.221M2. In the case of the Z2 spin liq- 
uid state we solve the full saddle point equations for k — 0.35 
and K = 0.5. Here k is a parameter measuring the strength 
of quantum fluctuations. In the large-A^ Sp(iV)theory the 



Owing to gapless fermionic spinons at nodal points, the 
U(l)-Dirac spin liquid state exhibits gapless spin-triplet exci- 
tations at the zero momentum point, F, and those momenta 
connecting two Dirac nodes, X. The lower bound of the 
two-spinon continuum has variations in energy approximately 
given by x J. Gapped bosonic spinons of the Z2 spin liquid 
state make the two-spinon spectrum also gapped with rela- 
tively large variations for the lower bound. On the other hand, 
the lowest spin-triplet excitation is completely non-dispersive 
in our valence bond solid theory. Since the flat dispersion of 
the lowest spin-triplet excitation is a distinct characteristic of 
the valence bond solid state, it may be used to distinguish our 
predictions from those of other spin liquid scenarios. 
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